Rationale for a new class of double-hybrid approximations in density-functional theory 
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We provide a rationale for a new class of double-hybrid approximations introduced by Bremond 
and Adamo [J. Chem. Phys. 135, 024106 (2011)] which combine an exchange-correlation density 
functional with Hartree-Fock exchange weighted by I and second-order M0ller-Plesset (MP2) corre- 
lation weighted by l^. We show that this double- hybrid model can be understood in the context of 
the density-scaled double-hybrid model proposed by Sharkas et al. [J. Chem. Phys. 134, 064113 
(2011)], as approximating the density-scaled correlation functional Ec[ni/i] by a linear function of 
I, interpolating between MP2 at / = and a density-functional approximation at Z = 1. Numerical 
results obtained with the Perdew-Burke-Ernzerhof density functional confirms the relevance of this 
double-hybrid model. 



The double-hybrid (DH) approximations introduced 
by Grimme [J, after some related earlier work 0, 
are increasingly popular for electronic-structure calcu- 
lations within density-functional theory. They consist 
in mixing Hartree-Fock (HF) exchange with a semilo- 
cal exchange density functional and second-order M0ller- 
Plesset (MP2) correlation with a semilocal correlation 
density functional: 



DH 



= a^Ef"" + {1 - a,)E.,[n] 
+{'i--ac)Ec[n]+a,E^^^, 



(1) 



where the first three terms are calculated in a usual self- 
consistent hybrid Kohn-Sham (KS) calculation, and the 
last perturbative term is evaluated with the previously 
obtained orbitals and added a posteriori. The two empir- 
ical parameters ax and ac can be determined by fitting to 
a thermochemistry database. For example, the B2-PLYP 
double-hybrid approximation [l| is obtained by choosing 
the Becke 88 (B) exchange functional Q] for Ex[n] and the 
Lee- Yang-Parr (LYP) correlation functional Q for Ec[n], 
and the empirical parameters a^ = 0.53 and ac = 0.27 
are optimized for the G2/97 subset of heats of formation. 
Another approach has been proposed in which the per- 
turbative contribution is evaluated with normal B3LYP 
orbitals rather than orbitals obtained with the weighted 
correlation density functional {1 ~ ac)Ec[n] 

Recently, Sharkas, Toulouse, and Savin have pro- 
vided a rigorous reformulation of the double-hybrid ap- 
proximations based on the adiabatic-connection formal- 
ism, leading to the density- scaled one-parameter douhle- 
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hybrid (DSIDH) approximation 
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{l^l)Ex[n] 
~l'Ec[n^/i]+fEr^ 



(2) 



where Ec[ni^i\ is the usual correlation energy functional 
evaluated at the scaled density ni//(r) = {l/l)^n{i/l). 
This reformulation shows that only one independent em- 
pirical parameter I is needed instead of the two parame- 
ters Ux and flc. The connection with the original double- 
hybrid approximations can be made by neglecting the 
density scaling 



Ec[ni/i] « Ec[n], 



(3) 



which leads to the one-parameter double-hybrid (IDH) 
approximation 



E. 



+ il-l^)Ec[n]+PE^^^. 



(4) 



Equation ^ exactly corresponds to the double-hybrid 
approximation of Eq. ([T]) with parameters Ox = I and 
Oc = P. 

Very recently, Brcmond and Adamo [T2j have proposed 
a new class of double-hybrid approximations where the 
correlation functional is weighted by {l~l^) and the MP 2 
correlation energy is weighted by P, instead of (1— Z^) and 
respectively. Applying this formula with the Perdew- 
Burke-Ernzerhof (PBE) [IJl exchange-correlation density 
functional, they have constructed the PBEO-DH double- 
hybrid approximation which performs reasonably well. 
In this work, we give a rationale for this class of double- 
hybrid approximations. For this, we start by recall- 
ing that the density-scaled correlation functional Edn^ /<] 
tends to the second-order Gorling-Levy (GL2) [ij] cor- 
relation energy when the density is squeezed up to the 
high-density limit (or weak-interaction limit) 
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FIG. 1: (Color online) Density-scaled correlation energy for the He (left) and Be (right) atoms as a function of /. 

Accurate calculations (from the parametrizations of Ref. are compared with different approximations: PBE without density 
scaling [Eq. Q], PBE with density scaling (from the parametrizations of Ref. |^), and linear interpolation between MP2 (with 
PBE orbitals) and PBE [Eq. ([9])]. The MP2 calculations for He and Be (including core excitations) have been performed with 
the cc-pV5Z and cc-pCV5Z basis sets [lol.[Tll|. respectively. 
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FIG. 2: (Color online) Mean absolute errors for the AEG (left) and BH6 (right) test sets as functions of the parameter I for the 
IDE, and DSIDH, and LSIDH approximations with the PBE exchange-correlation density functional. All calculations were 
carried out with the cc-pVQZ basis set. 



which is finite for nondegenerate KS systems. The 
GL2 correlation energy can be decomposed as (see, e.g., 
Ref.lil) 



£iGL2 _ ^MP2 



i^f"^ (6) 
where E^^"^ is the usual MP2 correlation energy cxprcs- 



E. 



MP2 



A 5Z $Z p 

ij ah 



\{(t)i<t>]\\(t)a<f)b)\ 



£b 



(7) 



with the antisymmetrized two-electron integrals 
{4>i4'j\\4'a4>b), a-nd E^^^ is an additional contribution 
involving the difference between the local multiplica- 
tive KS exchange potential and the nonlocal 
nonmultiplicative HF exchange potential 



E. 



AHF 



EE' 



(8) 



In both Eqs. ((T)) and (|8]), are the KS orbitals and 
Efc are their associated energies, and the indices i,j and 
a, 6 stand for occupied and virtual orbitals, respectively. 
The single-excitation contribution E^^^ vanishes for 
two-electron systems, and in most other cases is negli- 
gible so that the GL2 correlation energy is well ap- 
proximated by just the MP2 contribution (evaluated with 
KS orbitals), E^^'^ w E^^'^. This leads us to propose an 
approximation for Ec[ni/i] based on a linear interpola- 
tion formula 

E,[n,/i]^{l~l)Ef''' + lE,[n]- (9) 

Plugging Eq. (jH]) into Eq. we directly arrive at what 
wc call the linearly scaled one-parameter double-hybrid 
(LSIDH) approximation 

+{l-l^)Ec[n]+l^E^''^, (10) 

with the weights (1 — l^) and l^, thus giving a stronger 
rationale to the expression that Brcmond and Adamo 
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have proposed on the basis of different considerations. 
Further insight into this approximation can be gained by 
rewriting Eq. ([9]) in the alternative form 

EMi/i] « i?r' + 1 [EM] i?r') , (11) 

which can then be interpreted as a first-order expansion 
in / around I = with Edv] — E^^^ approximating 
the third-order correlation energy correction Ec \n] in 
Gorling-Levy perturbation theory. In comparison, the 
zeroth-order approximation i5c[«i/z] ~ E^^"^ plugged in 
Eq. ([2]) just gives the usual one-parameter hybrid (IH) 
approximation with the full correlation density func- 
tional [H, [13 

Elf' = lEf'' + il-l)EM+EAn]. (12) 

In this sense, the LSIDH approximation of Eq. (|T0| can 
be considered as a next-order approximation in I beyond 
the usual hybrid approximation. 

Figure [T] illustrates the different approximations to the 
density-scaled correlation energy i?c['^i/;] a-s a function 
of I for the He and Be atoms. The accurate reference 
curve is from the parametrization of Ref. |^. For I = 
1, it reduces to the exact correlation energy, while for 
Z = it is the GL2 correlation energy which tends to 
overestimate the correlation energy. In between these 
two limits, it is nearly linear with I. The PBE correlation 
energy without density scaling [Eq. ([3])] is the crudest 
approximation to Edrii/i]. The PBE correlation energy 
with density scaling (taken from the parametrization of 
Ref. gives a nearly linear curve. It is a fairly good 
approximation for He, but a less good approximation for 
Be where it underestimates the correlation energy for all 
I, especially at I = 0. This is due to the presence of 
static correlation in this system that is not described by 
the PBE functional. Finally, the linear interpolation of 
Eq. dH]) between MP2 (evaluated with PBE orbitals) and 
PBE appears as a good approximation for both He and 
Be. In fact, the linear interpolation is clearly the best 
approximation for Be, at least with the PBE functional. 

For a more comprehensive comparison of the different 
approximations, we have performed calculations on the 
AE6 and BH6 test sets ^ with the IDH, DSIDH, and 
LSIDH double hybrids for the PBE functional, using a 
development version of the MOLPRO 2010 program [19| . 
The AE6 set is a small representative benchmark set of 
six atomization energies consisting of SiH4, S2, SiO, C3H4 
(propync), C2H2O2 (glyoxal), and C4H8 (cyclobutane). 
The BH6 set is a small representative benchmark set of 
forward and reverse barrier heights of three reactions, OH 
+ CH4 ^ CH3 + H2O, H + OH ^ O + H2, and H -I- H2S 



~> HS + H2 . All the calculations are performed at the op- 
timized QCISD/MG3 geometries [Ig. We use the Dun- 
ning cc-pVQZ basis set [2l|, [22| ■ Core electrons are kept 
frozen in all our MP2 calculations. Spin-restricted calcu- 
lations are performed for all the closed-shell systems, and 
spin-unrestricted calculations for all the open-shell sys- 
tems. In Fig.[5J we plot the mean absolute errors (MAE) 
for the two sets as a function of I. For the AEG set, the 
MAEs of the DSIDH and LSIDH approximations are 
quite similar, and are both much smaller that the MAE 
of the IDH approximation for a wide range of I. For the 
BH6 set, the LSIDH approximation gives a MAE which 
is significantly smaller than those of both DSIDH and 
IDH for intermediate values of I. 

As regards the choice of the parameter I, the present 
data on the AH6 set gives an optimal value of I = 0.75 
for the LSIDH double hybrid, with a minimal MAE of 
3.59 kcal/mol, and for the BH6 set the optimal value is 
I — 0.70 giving a minimal MAE of 0.73 kcal/mol. How- 
ever, the MAE is not very sensitive to the value of / 
around the optimal value, and Bremond and Adamo [l3| 
have argued for using / = 0.5 in defining the PBEO-DH 
approximation, using a similar argument as the one used 
by Becke for his "half-and-half" hybrid j23j . 

In summary, we have shown that the new class 
of double-hybrid approximations named here LSIDH 
[Eq. PU)) ] can be understood as approximating the 
density-scaled correlation functional i?c["i//] by a linear 
function of I, interpolating between MP2 at ^ = and 
a density- functional approximation at Z = 1. Numer- 
ical results obtained with the PBE density functional 
confirms that the LSIDH approximation is a relevant 
double-hybrid model, and in fact tends to be more ac- 
curate than the DSIDH double-hybrid model [Eq. ©J 
in which density scaling is applied to the PBE correla- 
tion functional. More generally, it can be expected that 
the LSIDH double-hybrid model will be more accurate 
than the DSIDH double-hybrid model when applied with 
a density-functional approximation that is inaccurate in 
the high-density limit. We hope that this work will help 
constructing other theoretically justified double-hybrid 
approximations. 
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